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Abstract 



The classical limit of quantum q-oscillators suggests an interpretation of the 
deformation as a way to introduce non linearity. Guided by this idea, we 
considered q-flelds, the partition fumction, and compute a consequence on 
specific heat and second order correlation function of the q-oscillator which 
may serve for experimental checks for the non linearity. 



1. Introduction 

Recently the concept of quantum groups was developed in connection with 
the solution of nonlinear equations [1-5]. The special case of the quantum 
Heisenberg-Weyl group which produced the notion of quantum q-oscillator 
was studied by Biedenharn [6] and Macfarlane [7]. They introduced q- 
oscillators starting from generalized commutation relations containing an 
extra dimensionless parameter q. The physical sense of the quantum q- 
oscillator up to now was not clarified. The attempt to connect quantum q- 
oscillators with relativistic oscillator model has been discussed in [8]. The at- 
tempt to introduce q-oscillators and quantum group SU(2) in the frame of the 
generalized James- Cummings model has been discussed in [9]. Neverthless 
these attempts are based on pure mathematical properties of q-oscillators. 
The aim of our work is to show that the q-oscillator is a special case of a 
nonlinear quantum oscillator. 

The nature of q-oscillators of the electromagnetic field is clarified by the 
nonlinearity of field vibrations. The partition function for both classical and 
quantum q-oscillator is calculated. 

The spirit of the electromagnetic field q-quantum oscillator reminds us of 
such a model of electrodynamics as the nonlinear Born-Infeld model [10]. In 
addition to clarify the physical nature of quantum q-oscillator as the usual 
quantum oscillator with a "nonlinearity" of special type, the aim of this 
work is to study the classical q-oscillator. This turns out to be the classical 
nonlinear oscillator with the frequency depending on the amplitude by the 
particular functional dependence. We will show that in this sense the classical 
q-oscillator acquires a dependence on the amplitude of the physical pendu- 
lum. Assuming that the electromagnetic field be described by the q-oscillator 
(both -quantum and classical ones) we will study the physical consequences 
of this suggestion, namely how the parameters of the q-oscillators change 
the specific heat formula. This formula reflects the linear properties of the 
vibration of the vacuum at small amplitudes. It is clear that at very high 
intensities there must be influence of the nonlinearity of the vacuum vibra- 
tions. We will calculate the corrections to the partition function and to the 
specific heat formula. 

As for notations, we shall denote by h the deforming parameter, to avoid 
introducing a new symbol with respect to those already widely used in the 
literature [1-7]. We apologize with reader for this. We ask him to keep in 
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mind that Planck's constant is set equal to 1 everywhere in paper. 



2. Quantum q-oscillator 



Let us introduce the usual creation and annihilation oscillator operators a 
and at obeying bosonic commutation relations 

[a,<J] = l. (2.1) 

The classical dynamical variables to which a and a* correspond oscillate 
with a frequency uj — 1. It is known that the operators a, a\ 1 form the Lie 
algebra of Heisenberg-Weyl group. So, the linear harmonic oscillator may be 
connected with the generators of pure Heisenberg-Weyl Lie group. In view of 
the commutation relation (2.1) the usual scheme for generating the states of 
the harmonic oscillator is based on the properties of the Hermitean number 
operator n = a^a 

[a, n\ = a, [a\ n] 
Thus constructing the vacuum state |0) obeying the equation 



-at 



(2.2) 



a|0) = 



and the excited states 
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which are eigenstates of the number operator n 
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(2.3) 
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the matrix representation of the operators a and a) in the basis (2.4) have 
the known expressions 








VT 





... \ 










V2 
















V3 


V 








... J 


/ 











... \ 



















V2 







V 








... J 



(2.6) 



2 



while the number operator n is described by the matrix 

• \ 



n 



( 

10 

2 

v . . . 



(2.7) 



... J 



The Hamiltonian for such a system is defined as 

a* a + aaJ 
H = u u > 



(2.8) 



The q-oscillators may be introduced by generalizing the matrices (2.6) and 
(2.7) with the help of the q-integer numbers n q , 



n q = 



sinh rih 
sinh h ' 



q = e 



(2.9) 



Here h and q are dimensionless c-numbers, which appear at this purely math- 
ematical level. When H — 0, q=l and the q-integer n q coincides with n. 
Then, replacing the integers in (2.6) and (2.7) by q-integers we obtain ma- 
trices which define the annihilation and creation operators of the quantum 
q-oscillator, 



a n = 



( ^T q 
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a' = 
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... \ 
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(2.10) 



since the action of h q on eigenstates \n > is given by 

sinh nh , 



sinh h 



In > 



(2.11) 
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The above matrices obey the commutation relation 



[a q ,h] 



[a\,n] 



-J 



(2.12) 



but the commutation relations of the operators a q and a) q do not coincide 
with the boson commutation relations. Equation (2.1) is replaced by 



[a q , a J] = F(h) 
where the function F{n) has the form 

F{n) = f{h + l)-f{n) 



(2.13) 



(2.14) 



For h = (2.12) reduces to (2.1). 

In addition to the above commutation relation there exists the reordering 
relation 

a q a q — qa q a q = q~ n (2-15) 

which usually is taken as the definition of q-oscillators. 

It is worthy noting that the operators a q and a q can be expressed in terms 
of the operators a and a 1 " [11-13] 



where (see Eq. (2.9)) 



a f(n), aj = f{h)a\ 



(2.16) 

(2.17) 

(2.18) 

(2.19) 

In the Schrodinger representation the evolution operator of the harmonic 
oscillator 

(a^a + aaJ)t 



We have also 
and 



n a = 



a\a q 



[a q ,n q ] = F(n)a q 



H, fiq] = -a\ F (n) 



U(t) = exp 



—iuj 



(2.20) 



gives the possibility to find out explicitly linear integrals of motion which 
depend on time [14,15], 



e i0Jt a 



A(t) = U{t)aU- l {t) 
A\t) = U{t)(JU- l {t) = e- iu>t (J 



(2.21) 
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The matrices of the integrals of motion (2.21) in Fock basis may be obtained 
from the equations 



A(t)\n) = e ilJJt y/n\n-l),A\t)\n) = e~ i,JJt VnTT\n + 1) 

Let us now introduce the Hamiltonian 

a a at + a\a a 



H = u 



(2.22) 



(2.23) 



for which the evolution operator takes the form 
U q (t) = exp 



. (a q a\ + a\a q ) 
—lujt — - — - 



We have for the integrals of motion 

A q (t) = U q a q U q \ A* = U q a\U q l 
the following matrix expressions, 



(2.24) 



(2.25) 
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V ••• 

/ \ 

^ e -i{i q -%)^t o 

^Tge-^- 1 "^ 

V ••• / 



(2.26) 

These operators are the generalizations of the linear integrals of motion (2.21) 
to the case of nonlinear Hamiltonian. 

In general we can prove the following string of relations 



Ua a 



Ual 



a q U-\ 



(2.27) 
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where U± p is defined implicitely by expressing U as a function of the operator 

U± p (h) = U(n±p) (2.28) 

The integrals of motion A q (t), A q (t), which coincide with the operators a q , a q 
at t — 0, may be expressed in terms of the matrices (2.6) and (2.7) by means 
of the relations (2.26) and (2.27). 



3. Classical q-variables 

In this section we will use the inverse of the procedure proposed by Dirac to 
quantize a classical system, i.e. we consider a "dequantization procedure". 
We will obtain therefore a dynamical system on the classical phase space 
(q,p), or (a, a*), from the quantum system of the previous section, the latter 
then being the Dirac quantized version of the former. 

The usual harmonic oscillator vibrating with unit frequency may be de- 
scribed in terms of the following variables 



q + ip 
a = — ;=— 

q-ip 



a 



V2 ' 



with non-zero Poisson bracket 



and Hamiltonian function 



(3.1) 



{a, a*} = -i (3.2) 



H = aa* . (3.3) 
The flow associated with equations of motion is 

a = Ae^ t+ ^ a* = Ae~^ t+ ^ (3.4) 

where 

A = a(0)a(0)* 
2i a(0y 



1 i «(o) 

^ = o7 m 77n\7 (3-5) 
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For a frequency different from 1, the Hamiltonian is 

p 2 + u 2 q 2 

H = uj > (3.6) 

In terms of complex variables, we have now 

' ip i- \ 1 * f-iP , /- \ 1 



a = 



+ y/uq —= a* = —= + y/uq —j= 



H(a, a*) = ua a* 

{a, a*} = —i (3.7) 

From (2.16) we define the classical q-oscillator in terms of these new variables 
(a, a*), namely 



/ sinh Tiaa* / sinh Tiaa* 

Oi q = \ , - - a. a — \ , - - at (3.8) 

V aa* sinh n V olq* sinh n 

Their Poisson bracket can be computed and expressed in terms of themselves, 
obtaining 

{a q , a* q } = -*^^y / l + K| 4 (sinhft) 2 (3.9) 

so that we will consider a new system described by such q-variables with 
Hamiltonian function 

H(a q , a*) = ua q a*. (3.10) 
The equations of motion are then 

% 

smii n 

and complex conjugated, with solutions 



^1 + Icg^sinh/i) 2 ^ (3.11) 



a q (t) = a q (0)exp[-^^/l+ \ a q (0) | 4 (sinh^) 2 ] (3.12) 



-ituh 
sinh h 
and its conjugate. 

All this can be commented stating that we have performed a non canonical 
transformation, i.e. deformed the Poisson bracket structure, while preserv- 
ing the form of the Hamiltonian function and note that the above Poisson 
bracket has the standard bracket as a limit when H goes to 0. 
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Such a system can be, of course, re-written in terms of (a, a*) variables: in 
these coordinates the original Poisson bracket is unchanged while the Hamil- 
tonian function is 

sinh haa* 

H q (a,a) = uj — , - (3.13) 
sinh n 

This new dynamical system has a phase portrait which is the same as the 
usual linear harmonic oscillator. In facts, the equations of motion for the 
latter system, which evolves under the Hamiltonian H are 

a = —iua a* = iua* (3-14) 

and for the system with the Hamiltonian H q are 

a = —ioj q a a* = iu q a* (3.15) 

with 

— w — cosh/iaa* . (3.16) 
sinh n 

We notice that aa* is a constant of the motion for both systems. The main 
difference between such systems is that the frequency for the second one de- 
pends on the orbit while for the the first one is constant. 
This leads to interpreting q-oscillators as systems carrying a particular non- 
linearity. What in the harmonic oscillator is a constant u characterizing the 
evolution along any orbit, becomes here a function constant on each orbit 
separately, i.e. a constant of the motion for the undeformed harmonic oscil- 
lator. From the formulae (2.16) it follows a possible physical interpretation 
of the q-oscillator. The amplitude of the vibrations of the harmonic oscilla- 
tor is not connected with the phase of the vibrations which changes linearly 
in time. For the q-oscillator the frequency depends on its amplitude. The 
formula (2.16) reflects this classical nonlinear phenomenon for the quantum 
oscillator. Such a situation can be generalized to any linear system. 
We digress briefly to describe on what would be the situation in general. 
Let us have, indeed, a linear system 

X = AX (3.17) 

where X_ e R n and A a n x n matrix with constant entries. It can be 
integrated 

X{t) = e tA X(0) . (3.18) 
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Its integrals of the motion can be expressed in terms of the initial conditions 
X(0) = K . Replacing A with B(X) defined by 



where /j is any constant of the motion for the original linear system, the new 
system, which is non linear, 



can be integrated via exponentiation, because for any initial condition _B(X(0)) 
is a matrix constant along the trajectory originated from it. 
When the starting linear system is Hamiltonian, it turns out to be completely 
integrable. Therefore, eq. (3.20) is a generalization of our q-deformed oscil- 
lator. By employing various constants of the motion we could accommodate 
deformations with many parameters. There are other ways to get non linear 
systems out of linear ones, for instance by using a generalized reduction pro- 
cedure [16]. In the future we shall analyze and compare among them these 
procedures that use deformation or reduction to get non linear systems. 

4. Two dimensional q-oscillators 

In this section we will analyze a possible q-deformation of the two dimen- 
sional oscillator. There exists two essentially different possibilities to intro- 
duce the q-deformed creation and annihilation operators. They are based on 
our interpretation of the q-oscillators as describing nonlinear vibrations. If 
two modes in the linear limit have different frequences it is natural to deform 
the annihilation and creation operators for these two modes independently 
of each other, by simply using the formulae for one-dimensional q-oscillator 
discussed above. But if we want to take into account the isotropy of two 
dimensional oscillator and to conserve the U(2) dynamical symmetry even in 
the nonlinear vibration regime we might deform the two dimensional oscilla- 
tor in such a way that both constituent oscillators have the same frequency 
in the equal initial conditions for this regime too. Taking these remarks into 
account let us introduce two sets of operators 



(3.19) 



X = B(X)X 



(3.20) 



a q ± = a±f(n), aj± = f(h)a\± 



(4.1) 



9 



Here the number operator h is defined as follows 

n = n + + n_ (4.2) 

where 

n± = a±a± (4.3) 
and the operators obey the usual boson commutation relations 

[a±,<4] = l (4.4) 

all the others being zero. The commutation relations for the two components 
of the two-dimensional q-oscillator are then as follows 

[a q+ ,a q _] = 0, [a q+ ,a\_\ = a + al[f(n + 1) - f 2 (n)] (4.5) 

The other commutation relations are obtained from the conjugation of the 
above. The non-zero commutation relation of the annihilation and creation 
operators describing the same component is of the form 

[a q± , a\ ± ] = ala ± (f 2 (h + 1) - f(n)) + f 2 (n + 1). (4.6) 

Due to the nonlinear interaction of vibrations the operators a q+ and a q _ do 
not commute, as well as a q _ and a q+ . In the limit, when h goes to zero (q goes 
to 1) the written commutation relations lead to the usual boson commutation 
relations since in this limit the function / goes to l.The four operators 

a+a + ,a^a_,a+a_,a^_a + (4-7) 

form the generators of the U(2) group representation. This representation 
is realized in the space of states of two dimensional harmonic oscillator. We 
deformed the operator a + and a_ and they became the operators a q+ and a 3 _. 
But the deformation given by the formulae (4.1) is such that the operators 

al + a q+ , al_a q -, a q _a q+ , a\_a q+ (4.8) 

realize the same representation of the same undeformed group U(2) in the 
space of the same two dimensional harmonic oscillator states. This happens 
because these new generators differ from the former ones only by the coeffi- 
cients which depend on the Casimir operator of the rotation group. So, for 
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any irreducible subspace these coefficients are simply c-numbers and the op- 
erators (4.10) close on the Lie algebra of the group U(2) in these subspaces. 
It can be shown that the following re-ordering relation holds 

a q+ a\ + + a g _<4_ - q{a\ + a q+ + a\_a q J) = q~ n (4.9) 

Had we q-deformed the two components of the oscillator as independent ones, 
namely as those appearing in the previous section 2, the right hand side of 
the Eq.(4.11) would be replaced by the term q~ n + + q~ n - . This is a different 
instance of the non linearity of the q-oscillators. 

If we deal with many more oscillators, we can use mixtures of these q- 
deformation procedures according to which linear transformations we want to 
survive after the deformation process. Let us clarify this statement. We now 
study transformations between two-dimensional q-oscillators which are anal- 
ogous to canonical transformations of linear oscillators, (a±, a±) and (b±, b±), 
by requiring that 

a + a\ + a-al = b + b{ + b-bl . (4.10) 
The resulting linear transformations are 

a± = A±b + + B±b- (4.11) 

4 = A* ± b\ + B* ± bl (4.12) 
with A± and B± complex numbers such that the matrix 

A + B, 
A_ B_ 

belongs to U(2). Substituting the above formulae in (4.1) and (4.2) we have 
consequently 

a 9 ± = (A±b + + B ± b-)f(n) (4.13) 
al ± = f(n)(A* ± bl + B* ± bl) (4.14) 
Here the number operator n is written as follows 

h = n b+ + n b _ (4.15) 

where 

n b+ = b{b + (4.16) 
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n 6 _ = bib- . (4.17) 
If we consider a dynamical system described by the Hamiltonian, 

H = ^( a l+ a q+ + 4- Q 9- + a q+°i+ + a q- a l-) ■ ( 4 -18) 

Written in terms of the number operator, H acquires the following form 

1 



H =2 



hf(n-l) + (h + 2)f 2 (h + l) (4.19) 



which is manifestly invariant under the previous transformations. In the Fock 
basis | n+ > |ra_ > the eigenvalues of H are immediately computed as n has 
positive integers eigenvalues and each of them has degeneracy n + 1. 
It is obvious how to define the variables for the classical 2-dimensional q- 
oscillators. For instance 

a q ± = a±f(n) (4.20) 
where (a±,a±) are variables for the usual harmonic oscillators and 

n = n + + n_ (4.21) 

with 

n± =| a± | 2 . (4.22) 

The computation of the non-zero Poisson brackets then gives the following 
result 

(hn) cosh nh — sinh rih 
n 2 sinh h 

and 



{a q ±, a* qT } = -ia±a% — (4.23) 



{a„±, «*+} = r~r[(l ) sinhn/i + hn± cosYinh] . (4.24) 

q nsinhn, n 

Thus, it is clear that those transformations which preserve the term respon- 
sible for the non linearity will be compatible with the q-deformation. More 
explicitely, this means that these transformations commute with the defor- 
mation procedure. 



5. q- fields 
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We take into consideration the usual expansions of classical and quantum 
fields on an interval of length L = 1 in terms of usual harmonic oscillators. 
For simplicity we start considering a classical real massless scalar field <f>(x, t) 
and its conjugate ir(x,t), satisfying also <p = n 

— oo,+oo — oo,+oo 

<Kx,t)= £ e ik ^ qj (t) ir(x,t) = £ e-*>* Pj (t) . (5.1) 

3 3 

Notice that here qj and P j are complex variables which anyhow can be ex- 
pressed in terms of harmonic oscillators as, for instance in [17] so that also 
in this case the usual Poisson bracket holds 

{qi,Pj} = $ij kj = 2nj. (5.2) 

The use of such fields is aimed at a study of some aspects of electrodynamics 
with massless q-photons. Following [17], we have 

j V ZUJ j 
n(x,t) = -ij:e-^M(a_ 3 -a*) 

3 V ^ 

U3j = \2TTj\ . (5.3) 

These variables (<x,-, a*) with j — ± integers can be written in terms of usual 
oscillator variables for positive and negative frequencies (aj±, a*±) with j 
positive integers, namely 

% = . (5-4) 

a .j = (5.5) 

and their complex conjugates. The Poisson brackets of the usual oscillators 
are of the form 

{a j± ,a j/± } = , (5.6) 
{a j± ,a*, T } = , (5.7) 
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and 

{a j± ,a* jl± } = -iSjj. . (5.8) 

We have 

{(f)(x,t), ir(x',t)} = 5{x-x') . (5.9) 

We are now in a position to define a pair of classical q-fields, noting that 
there is not a unique definition due to the fact that the frequency depends 
on aja*. 

By using previous equations, we write 



Tr q (x,t) = -iJ2e- i2 ^V^ jq 



(5.10) 



V2 



where 



"» = IMn Un ' < 5 ' n > 

and aj q and a* q are given as functions of aj± and its complex conjugate; for 
instance, 



sinh( \aj + 


2 + 


a 3~ 


2 )^ 


(K+ 


2 + 




2 ) sinh/i 



(5.12) 



The computation of the Poisson brackets for these fields gives results which 
are different from the usual ones, as it was expected, up to corrections of 
the order of h 2 . The above fields have been constructed allowing a coupling 
between positive and negative frequency modes belonging to the same index 

j- 

To give an idea of the effects of a q-deformation, we give explicitely the 
behaviour of the Poisson bracket with H — > between conjugated fields, 
which have been q-deformed in another way that allows for much simpler 
computations. Such new fields are defined as follows 

U*> t) = E e^(^= + -§M (5.13) 
ir q (x, = -< E e-^i^a^ - ^fa* q ) (5.14) 
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where 

h 

Ujq = '^'imhli cosh (W 2/l ) ( 5 - 15 ) 
and the q-deformed oscillators are each other independent and defined as in 
(3.8). That is we are allowing only the self coupling. We can compute the 
Poisson brackets between such fields at equal times and obtain 

h 2 

{<f> q (x,t),n q (x',t)} = 5{x - x')(l - — ) + h\ Cl + c 2 ) (5.16) 



where 



and 



c i = 7 EK^'^d^f + + cc] (5.17) 



i>o 



c 2 = ^E[(^ W) («?l«il 4 -«-il«-il 4 )) + c.c.] • (5.18) 
It can be shown that 

{0 9 (x,t),0,(x',t)} = O(/i 2 ) (5.19) 

{n q (x,t),Tr q (x',t)} = 0(ti 2 ) . (5.20) 

All this gives an account of the effect of q-nonlinearity as far as the principle 
of locality is concerned. The corrections destroy the ^-function behaviour of 
the Poisson bracket between conjugated fields and are proportional to the 
fourth power of the amplitudes of nolinear modes vibrations. We shall come 
back to these aspects in the future. The quantized fields are defined from eqs. 
(5.13-14) replacing c-numbers (aj q ,a* q ) with the creation and annihilation 

operators (a,j q , aj ? ) with ujj q becoming also an operator LUj q , namely 

fax, f) = £ e^(-J=a, 9 + -T^=al 3q ) (5.21) 

3 \/ 2uJ j q \/2V-jq 

and 



with 



ir q (x, t) = £ e«( f^a_ iq - ffa) q ) (5.22) 



3 9 1 h) 
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Such a definition of quantized q-fields just follows the usual prescript- ion, 
particularly for what concerns the j = term in the sum and holds since the 
operators 0Uj q are invertible for j ^ .0. Their commutation relations at equal 
time can be calculated and as it was easily expected, we have that such fields 
are not local. 

It should be noticed that we have not addressed the problem of deformation 
of the differential Maxwell equations, i.e. we are not exhibiting a q-deformed 
d'Alembertian. As a matter of fact, at this stage we do not even know 
whether our q-fields are solutions of some differential equation. These and 
other related aspects will be taken up in a forthcoming paper. 

6. Partition and time-corrrelation functions 

The partition function for a single q-oscillator 

oo 

Q(/3,70 = £exp(-/3n,) (6.1) 

n=0 

is a function of the temperature, through the variable (5 = u/kT, and the 
parameter h. In particular, % is very small and the series representing Q 
can be represented in analytic forms. We will consider first the limiting 
expression of Q obtained by letting % — > for fixed (5 thus obtaining a series 
representation of the form, 

Q(P,h) = Q + 1 -h 2 ^Q + ... (6.2) 
2 an 

where Qo(P) = Q{P,fy = 1/(1 — e_/3 ) represents the partition function of 
the ideal harmonic oscillator. In particular, we will see that the above series 
does not converge for (5 < h. In view of this, we will discuss the limit of Q 
for h — > for given value of the ratio ( = {3/ sinh h. The quantum partition 
function obtained in this second case is well approximated by the classical 
one. 
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6.0.1 Partition function for given f3 



When we keep (5 fixed and let h go to zero we can represent Q as a power 
series in Ti by relying on the power expansion of n q , 



= n + — (n d - n) + ^- [n 5 - yn d + -n) + • • • (6.3) 
Then, 

oo 

Q(P,h) = ex P ( ~P n q + P n ) ex P ( ~l3n) 

n=0 

( w 2 \ 

« Qo(/3) ( 1 + ^-/9 3 (Q 3 - Qi)) 

~ W) (i + ^ 3 (g 3 -Qi) 

+ W A (^ 6 (Q 6 - 2Q 4 + Q 2 ) + ^ 5 (Q 5 - yQ 3 + ^Qi; 
where W — fi/f3 and 

Qo(0) 



(6.4) 



1-e-/ 3 



Qi = i-Q 

Q2 = l-3Q + 2Ql 

Q 3 = 1 - 7Q + 12QI - 6Q 3 (6.5) 



Accordingly, we have 

Q(P, h) » Q (/5)(l - W 2 /? 3 Qo(Qo - I) 2 ) (6.6) 
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Consequently, the average energy E and the specific heat C are respectively 
given by 

E = -^ = E + h 2 ^([3Q (Q -l) 2 ) 

C = C -(3 2 h 2 ^- 2 (PQ (Qo-l) 2 ) (6.7) 

It is worth discussing the behavior of the above expression of Q for given 
W and j3 — ► 0. Since for small (3 Q m 1/ [3, then 

Q(P) « Q (l-W 2 + W 4 + ...) 
E{(3) « £ -2W 2 

C(/3) « C -6W 2 (6.8) 

Accordingly, the above representation of Q is valid only for (3 2 > 6h 2 , that is 
for W < 1/y/E. 



6.0.2 Partition function for given f3/% 

For given ( = (3/U we can rely on the following inequalities for the partition 
function 

/ e-^ {x) dx < Q(f3, h)<l + / e~ pE " (x) dx (6.9) 
jo Jo 

where E q (x) is a monotone function of x which coincides with n q for x = 
n. The above integral, which represents the classical limit for the partition 
function, can be represented as a combination of Weber function E and 
Neumann function N [18], 

fCO 77- 

I e-MMdx = ~[N (0 + Eo(0\ (6.10) 
being ( = (5/ sinh H. In particular, 

E(0- f (~ 1 )"(C/2) 2 " +1 f6U x 
^ (0 "~^ ^(n + 3/2) (6 - U) 
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It is noteworthy that for fixed ( there exists a limiting value h of h such 
that for h < h the partition function is so large to coincide with the classical 
one. In fact, for ( very small we have 

-W^ + N ^ X 2r4w( +0(C2> (6 ' 12 ' 

Hence, for ( < .1 the classical and quantum partition functions are almost 
coincident for ( > 10h. 

In conclusion, while the quantum partition function of the q-oscillator is 
well approximated by the classical limit for h going to zero, this is not true 
for the classical harmonic oscillator. 

Next, we consider the average energy E for which we can write the fol- 
lowing inequalities 

E < <E<E > (6.13) 
where the lower (-E<) and upper (E>) bounds are given by 

Plugging (6.10) into the above expressions yields 

(2smhhW)/(ne)-E' (()-N' (() (2 sinh hW)/(ne) - E' (£) - N' (() 

-Ah/n + E (C) + N (C) E (C) + N (C) 

(6.15) 

When % is so small to satisfy the following inequalities 

— ^EoiW' 1 ) +N (W' r ) 

7T 

« _ N ' o{W -i) (6.16) 

7re 

we can approximate E by 

E ^_^K±K (617) 

sinh h E + N v 1 
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Next, letting (3 — > into the above equation we have 



(3E&-—— (6.18) 



In 

Consequently, the specific heat decreases for T — > oo as 

CK hTr < 6 ' 19) 

Thus the behavior of the specific heat of the q-oscillator is different from 
the behavior of the usual oscillator in the high temperature limit. This 
property may serve for an experimental check of the existence of vibrational 
non linearity of the q-oscillator fields. 



6.0.3 Time-correlation function 



For characterizing adequately the thermalized q-oscillator it is worth dis- 
cussing the relative second-order correlation function j q (P, h; t). Defining 



Q(P, M) = E exp(-Pf 2 (n) + itF(n - 1)), 

n=0 



we have 



lq (z, h]t)=Tr(paUt)a q (0)) = 



1 dQ{P,h;t) 



(6.20) 



(6.21) 



To compute j q (P, h; t) for h — > we can use the approximate expression 



of Q(P,h;t) 
Q(P,h;t)ne u 

where 



^ ,n\ +2 n< n >-<n> . < n > - < n z > 
Q (P) - h 2 (3 + it 



< n k >= J2 n k e~ pn A; = 0,1,.... 



(6.22) 
(6.23) 



n=0 



It is worthwhile noting that 7 g depends on the hierarchy of photon dis- 
tribution momenta < n >,< n 2 >,.... This circumstance marks the main 
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difference between q-oscillators and standard ones. In fact, the dependence of 
the correlation function on the intensity is a general property of q-oscillators, 
holding true also when the q-oscillator is in a generalized coherent state. This 
property could be used for testing experimentally the possibility that elec- 
tromagnetic fields behave like q-oscillators. In fact, accurate measurements 
of the dependence of the field correlations (second and higher orders) on the 
field intensity could help in finding upper bounds for the small quantity h. 
For t — the correlation function gives the mean value of h q ,i.e. 

n q = lq ((3,h;0). (6.24) 

Such a distribution function is expressed in terms of Eq. (6.21) in which 
Q(P, h] t) is taken at t=0. 

q-deformed Bose distribution n can be obtained by the same method 
starting from the Hamiltonian operator (2.23) with u — 1 and not from 
H = a q a q and one obtains 

™ = ^o-y ^((n 2 )o- (n)l) + ^((n 3 )o-n (n 2 ) ) + ^((n i ) -n (n 3 ) ) , 

(6.25) 

in which n is the usual Bose distribution function and 

n oo 

n k = 2 sinh '- ]T n * e -PW>) . (6.26) 

2 n=0 



7. Conclusions 

We conclude this paper by pointing out some ideas of the work. 
As we understand now the one-dimensional q-oscillator is nothing else than a 
nonlinear oscillator with a very specific type of the nonlinearity. Namelyits 
frequency depends on its energy as hyperbolic cosine of the energy. Thus, 
classical motion of such nonlinear oscillator is described as motion of q- 
oscillator. So, in this case the frequency of the vibrations (or the velocity 
of motion) increases exponentially with the energy. After standard quanti- 
zation we obtain quantum q-oscillator which is nothing else than a nonlin- 
ear quantum oscillator with specific anharmonicity described by the infinite 
power series in energy. From the discussion above it is now clear that there 
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are different approaches to generalyse this picture. In the frame of nonlinear- 
ities, even when the frequency of vibrations depends only on the energy we 
could choose other functional dependencies introducing functions different 
from hyperbolic cosine. We could call these oscillators f-oscillators where 
/ is now the function determining the dependence of the frequency on the 
energy.This function is related to the Poisson brackets of new a /-coordinates 
where af = otf(aa*) and the a is the harmonic oscillator coordinate. Of 
course, this function has to contain a dependence on the parameter H such 
that for H — > the function / — > 1. 

The quantum /-oscillator obtained by canonical quantization of the above 
is again a nonlinear oscillator with a different dependence of u> on the en- 
ergy. We could also introduce other nonlinearities ny making the frequency 
to depend on other constants of the motion, different from the energy. This 
gives the possibility to undertake a classification of such nonlinearities and 
especially it is worthy for multidimensional oscillator. Here we may introduce 
a deformed nonlinear motion considering the frequencies of the oscillators de- 
pending either on their own energies or on energies of some groups of other 
oscillators. So, to preserve the U(N)-symmetry of multidimensional isotropic 
oscillator we may introduce the frequency dependence on energy to be a 
function of total energy of all the oscillators. If each oscillator is deformed 
as independent one we destroy the U(N)-symmetry of the initially isotropic 
multidimensional oscillator. 

In the quest for adequate field theories, we limited ourselves to consider 
some possible effects in systems which can be represented as a collection of 
single oscillators. In particular, we showed how the specific heat formula 
is changed by the q-non-linearity. In addition, the second-order correlation 
function depends on the average occupation number of the q-oscillator. This 
property marks the main difference with the standard oscillator. 

In spite of mathematical beauty of q-deformation procedure, (and /- 
deformation too) the most important issue is to envisage possible experi- 
mental tests of the possibility of describing electrodynamic phenomena by 
means of these q-oscillators. At a first thought , non-linear optics exper- 
iments appear as potential candidates for obtaining upper bounds for the 
value of h. However, these experiments require the development of a q-field 
formulation of QED. For example, if a laser beam behaves like a collection 
of q-oscillators interacting with the matter in the standard way, an opti- 
cal rectification experiment carried out at different laser intensities could be 
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used for testing the deviation of the boson spectrum from the ideal equis- 
paced pattern. As an instance, apparatus presently developed for carrying 
out accurate interferometric experiments (f.i. LIGO in USA and VIRGO in 
Italy-France) on gravitational waves, in view of their ability to reducing the 
photon noise limit by averaging over periods of the order of one year could 
in future become ideal candidates for performing these tests. 

Another class of experiments could be based on the dependence of time- 
correlation function of q-oscillator on the intensity. The interpretation of 
these experiments does not depend critically on the development of a q- 
field QED. In fact, these measurements can be carried out in vacuum. For 
example, a continuous laser beam could be split in two beams of different 
intensities and sent to two spectra analyzers. Measuring the widths of the 
two spectra it would be possible to put an upper limit to the deformation 
parameter h. These measurements could last for years thus guaranteing 
a signal-to-noise ratio adequate for appreciating very small values of the 
parameter. 

The non-linearity may also produce sqeezing [19] or may be related to 
non locality and lattice structure in space and time. The q-oscillator and 
squeezing has been discussed in [20]. The non locality in time has been 
studied in [21] where a q-deformed linear version of Klein-Gordon equation 
has been suggested. We would like to point at the recent work [22] where 
similar ideas are discussed. The connection of q-oscillator with anharmonic 
oscillator has been also discussed in [23]. We will present a more detailed 
analysis of possible generalizations and experimental consequences in other 
future publications. 
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